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Abstract 
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quark pair production diagrams with subsequent top quark decay into a bottom quark and a charged Higgs 
boson. We compare the NLO differential cross sections obtained in the bottom partem picture with those 
for the gluon-initiated production process and find good agreement. The effects of supersymmetric loop 
contributions are explored. Only the corrections to the Yukawa coupling are sizable in the potential dis- 
covery region at the CERN Large Hadron Collider (LHC). All expressions and numerical results are fully 
differential, permitting selections on the momenta of both the top quark and the charged Higgs boson. 
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I. INTRODUCTION 



The elucidation of electroweak symmetry breaking is an important goal of particle physics. In 
the standard model, one neutral scalar Higgs boson is assumed to exist, and it is associated with 
the generation of the masses of the electroweak gauge bosons and of the fermions. The neutral 
Higgs boson has not yet been observed, and direct searches place a lower limit of about 114 GeV 
on its mass [1]. Extensions of the standard model include the possibility of more Higgs fields. 
The minimal supresymmetric standard model (MSSM) requires two doublets to give mass to up- 
type and down-type fermions and to cancel anomalies. The doublets yield five physical Higgs 
bosons: two neutral CP-even states, a CP-odd state, and a pair of charged scalars. At lowest 
order in perturbation theory, the masses and couplings of these states depend on two parameters 
which may be chosen as the pseudoscalar mass rriA and the ratio of the two vacuum-expectation 
values tan/3 = v 2 /v 1 . Comprehensive analyses have been performed of the expected coverage 
of the (tan/3, m^) parameter space at the CERN Large Hadron Collider (LHC) J2J]. While the 
observation of at least one of the two CP even Higgs bosons may not pose a problem for the 
LHC 1 3], it will be challenging to distinguish it from its standard model counterpart over a large 
fraction of the parameter space. For small values of tan (3 the only viable channel in which to 
observe a heavy Higgs boson could be the resonant production of the scalar H with subsequent 
decay to hh — > 6677 |4|], where b is a bottom quark. For large values of tan (3, the identification 
of a charged Higgs boson would provide evidence for a Higgs sector beyond the standard model, 
meaning at least two Higgs doublets, and possibly a supersymmetric Higgs sector. 

If the charged Higgs boson is lighter than the top quark t, there is a good chance that it will be 
discovered via the decay channel t — > bH + in pp collisions at the Tevatron collider at 2 TeV, or in 
pp collisions at the LHC at 14 TeV. Searches in the Run I data samples by the CDF and DO col- 
laborations at the Tevatron [Q] place significant bounds on the mass m H and tan (3. If the charged 
Higgs boson is heavier than the top quark, then its observation at hadron colliders becomes more 
problematic. In particular, there is no tree level coupling of a single charged Higgs boson to gauge 
boson pairs, and the production of H ± is inaccessible in weak boson fusion. The cross section 
for H + H pair production is likely to be too small, and the heavy quark backgrounds may be too 
large for the observation of charged Higgs boson pairs, unless additional supersymmetric parti- 
cles enhance this loop-induced rate [Q, Q]. The situation is similar for the associated production 
of a charged Higgs boson with a W boson. The standard model leads to a fairly small rate, but 
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supersymmetric particle loops might enhance the rate considerably 

The most promising search channel for a heavy H ± is the associated production of a top quark 
and the charged Higgs boson pp — > tH~ + X and pp — > iH + + X via the intermediary of a bottom 



quark coupling [ 10]. Throughout this paper we present results only for the tH~ channel, unless 
stated otherwise. If both the tH~ and iH + channels are included the rates increase by a factor of 
two. Advanced detector simulation studies have been done for the decay channels H — > tb Jill 
and H — » rv Jl^l . The advantage of this production mode is that the Yukawa coupling to a top 
quark and a bottom quark is enhanced by a power of tan (3 for large values of tan (3, as are the 
bottom quark Yukawa couplings of the heavy neutral Higgs bosons. 
We consider the leading order partonic subprocess to be 

gb - tH~, (1) 

with the initial state bottom quark taken as a constituent of an incident proton. The set of next-to- 
leading order subprocesses includes partonic reactions such as gg — > btH~ . The total rate for the 
process pp — > btH~ + X receives large corrections from collinear logarithms, originating from 
the radiation of a forward bottom quark jet I13I [ul . These logarithmic terms can be resummed 
to all orders in the strong coupling strength a s , leading to the bottom parton picture |13] with an 
appropriate bottom quark factorization scale 11141 This resummation of large collinear 

logarithms is valid not only for charged Higgs boson production, but it is generic as long as there 
is a large mass scale M that provides \og(M /pr,b) behavior. In our case M = m t + rriH- The 
comparison of higher order predictions for total cross sections of neutral Higgs boson production 
shows impressive agreement between gluon-initiated and bottom parton results Jnl [lsl . but the 
absence of a heavy scale in neutral Higgs boson production in association with bottom quarks may 
have an impact on some final state distributions. 

The tH~ production cross section can be evaluated with or without integration over the phase 
space of the final state bottom quark, corresponding to whether an accompanying final state bottom 
quark is observed or ignored. The term "exclusive" is generally used to refer to a situation in which 
the final state bottom quark is observed, and "inclusive" is used to refer to the case in which the 
final state bottom quark is ignored. 

Calculations of the next-to-leading order (NLO) total cross section for pp — > tH~ + X are 
available in perturbative quantum chromodynamics (QCD) as well as in supersymmetric 

QCD J3,0]. If we use the same choice of parameters, most importantly the same renormalization 



and factorization scales, and the same scheme for renormalization for the Yukawa couplings, our 
numerical results are in good agreement with those of Ref. [19]. We comment in greater detail 
in Sec. HU C about the differences in scheme choice. The inclusion of NLO contributions to the 
gb initial state process merges the "inclusive-type" process gb — > tH~ with the "exclusive-type" 
process gg — > btH~, whose contribution appears as part of the set of NLO diagrams. The NLO 
contributions increase the reliability of the theoretical predictions by reducing the renormalization 
and factorization scale dependence of the total rate. 

In this paper, we present fully differential NLO cross sections for the process gb — > tH~. The 
differential distributions are desirable, as are predictions of expected correlations among the final 
state observables, since selections on final state kinematic variables must be made in experimental 
studies, for reasons of event acceptance and background rejection. In Sec. HD we outline the two- 
cutoff phase-space slicing method which we adopt. We then study the NLO production rates at the 
LHC and the Tevatron with the associated theoretical uncertainties. We present typical kinematic 
distributions and momentum correlations in Sec.|ni| For searches in the framework of the MSSM 
we examine the effects of leading and sub-leading supersymmetric QCD corrections in Sec. IIVI 
Conclusions are summarized in Sec.fVl 



The results in this paper go beyond those of Refs. II14LI19I1 in several respects. The calculation 
presented in Ref. 111411 uses a one-cutoff method for NLO calculations. The two-cutoff method 
in this paper permits a fully differential treatment of the final state particles. In this paper, we 
show that the bottom parton approach is justified for differential cross sections as well as for total 
cross sections. A more extensive discussion of cutoff dependences is presented in this paper in- 
cluding two-dimensional plots. In our treatment of SUSY-QCD corrections in Sec.|lVj we include 
an evaluation for the Snowmass points and slopes (SPS) parameters and the impact of resumma- 



tion of the corrections. Going beyond Ref. 11911 . we include SUSY-QCD corrections, and an 



exploration of the interplay of bottom parton and gluonic contributions. The matching of cross 
sections for small charged Higgs boson masses is new. Matching near the top decay threshold is 
of considerable interest for LHC experiments. 

II. NEXT-TO-LEADING ORDER QCD CORRECTIONS 

Throughout the paper we use a running strong coupling a s , a bottom quark Yukawa coupling 
Ub, and a top quark Yukawa coupling y t consistent with the order of perturbation theory, i.e. one- 
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Figure 1 : The leading order Feynman diagrams for the production process gb — > tH . We indicate how the 
bottom partons are created through gluon splitting. 

loop running for the leading order and two-loop running for the next-to-leading order results. 
If not stated explicitly otherwise we neglect the bottom quark mass m& in the phase space as 
well as in the matrix elements, while naturally keeping it in the bottom quark Yukawa coupling. 
The Yukawa coupling is normalized to mft(m&) = 4.2 GeV and a bottom quark pole mass of 
4.6 GeV. Moreover, we use the CTEQ5 parton densities ll^lll . We refer to the K factor defined 
by K = onlo/olo* me rat i° °f me NLO cross section over the leading order cross section. The 
default choices of the renormalization scale and the factorization scale are taken to be proportional 
to the hard scale in the process 

fx° R = M/2 n° F = M/5 (M = m H + m t ). (2) 

A detailed explanation of these choices may be found in Sec. HUB. 

At leading order in QCD we start from the parton-level production process in Eq. <[T]), with 
the diagrams depicted in Fig. [T] It is appropriate to define the bottom quark as a parton at high 
energies since the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation [22] 
resums large logarithmic log pr,b contributions from small p™™ ~ m& to a maximum value p™^ x 
(which in turn determines the 6-quark factorization scale). The bottom parton density is therefore 
not suppressed by a simple power of a s . The leading order cross section for a process involving 
one incoming bottom parton and an incoming gluon is of order a s y% t , where y^ t is the bottom- 
top-Higgs Yukawa coupling. If y^ t is written as y\ t = yf + y\, y t and y b are terms proportional 
to (m t / tan (3) and [m b tan (3) respectively. For large values of tan (3, y^t is dominated by y b . The 
validity of our choice of the process in Eq. CQ) as the leading contribution is confirmed by our 
numerical results, namely that the perturbative series is well behaved over a wide range of scales. 

There are two classes of NLO contributions: 

(1) Virtual gluon exchange corrections to the lowest order process and the corresponding real 
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gluon emission corrections, both of order a 2 s yl t , 

gb — > tH~ (virtual correction) 

gb - tH-g. (3) 

(2) The purely gluon-initiated and the purely quark-initiated diagrams, which lead to cross sec- 
tions also of the order a 2 s y 2 t , 

gg, qq, bb — > tH~b, bb — > tH~b, 

bq tH'q, bq -> tH~q. (4) 

Because we neglect the bottom quark mass in the phase space and in the matrix elements, the 
purely gluon and purely quark initiated subprocesses are divergent in the collinear limit. In our 
calculation these divergences are removed through mass factorization, i.e. the proper definition of 
all parton densities at NLO. 

One may think about an alternative treatment of the associated production process, namely 
to start with the process gg — > tH~b as the leading contribution. These diagrams are part of 
the a s correction to the bottom-gluon fusion process. For a choice of the factorization scales 
^F, g = ^F,b — ► m-b the bottom parton density vanishes, in contrast to the gluon density, which is 
stable and well defined down to scales of the order of Aqcd- In a physical picture of this limit we 
consistently switch off all large collinear logarithmic contributions, because p™^ x = iip,b- We are 
then left with only the purely light-flavor qq and gg induced processes listed in Eq. @. We use 
this limit of a small bottom quark factorization scale in Sec.|ni|to check the impact of the bottom 
parton picture on the final-state differential cross sections. 

A. Phase-space slicing 

One of the main tasks of our calculation is to integrate the three-body matrix elements over 
the phase space of the unobserved particle in the final state. The situation is different from the 
case of the single particle inclusive calculation in which one integrates over the phase space of two 
particles in the final state. We wish to retain control over the kinematic variables of a second par- 
ticle in the final state, while at the same time integrating over enough of the phase space to ensure 
cancellation of all infrared and collinear divergences. Several techniques have been introduced for 
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these purposes. The phase-space slicing J23IQ] and the subtraction methods [ 25 ] are two ways to 
extract the singularities in the real emission contributions as exclusively as possible. All relevant 
information needed to compute a 2 — ► 2 particle NLO cross section with the two-cutoff slicing 
method is compiled in Ref. I24L We follow this description closely in our calculations. 

The ultraviolet divergences in the virtual 2^2 corrections are handled with dimensional 
regularization. The heavy final state masses are renormalized in the on-shell scheme, while all 
couplings — the strong coupling as well as the bottom quark and the top quark Yukawa couplings 
— are renormalized in the MS scheme. The mismatch of the pre-factors between the virtual 
corrections and the counter terms leads to the usual explicit log fiR dependence of the NLO cross 
section on the renormalization scale. The choice of the renormalization scale and that of the 
factorization scale jj F are discussed in Sec. Ill Bl 

Virtual gluon exchange and real parton emission lead to both soft and collinear divergences. 
They are extracted with dimensional regularization and partially canceled with each other and 
partially removed through mass factorization, i.e. the consistent definition of parton densities. The 
situation is relatively simple for processes that have different initial states from the leading order 
gb case because no soft divergence appears. Schematically, we can write the contributions arising 
from the processes in Eq. @ as: 

= d*?% + do*% + daf MF . (5) 

The label HC indicates hard collinear divergences, which cancel with the universal contributions 
from hard mass factorization (HMF). The collinear phase space region, which appears for all 
2^3 kinematics, is defined as the region in which the value of the corresponding invariant for 
the two possibly collinear momenta p, and pj falls below (pi+Pj) 2 < 5 C s, where y/s is the partonic 
center of mass energy. The squared matrix element is finite in the non-collinear phase space region 
and can be integrated numerically, creating an implicit logarithmic dependence on the cutoff 5 C . 
Applying the hard mass factorization corrections, we subtract the hard-collinear contributions as 
a convolution of the leading order matrix element and the appropriate finite splitting function, 



multiplied by log5 c [24]. The mismatch of pre-factors leads to an additional explicit dependence 
of the NLO cross section on the factorization scale log fip. 

The situation for virtual and real gluon emission in the gb initial state is slightly more involved 
because additional divergences appear due to soft gluon emission. Soft gluon emission is defined 
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by the non-invariant gluon energy constraint E g < S s y/s/2. The cross section can be written as 

da g = da s virt + doS& + da s 2 ^ 3 + da s 2 Zs + da SMF + daf% + da^ g + daf MF , (6) 

where the label S means soft, SC soft-collinear, and SMF soft-mass-factorization. The additional 
soft 1/e and overlapping soft-collinear 1/e 2 divergences appear in the virtual corrections as well 
as in the real gluon emission corrections. The divergences cancel among the virtual correction, 
the real gluon emission correction, and the contributions from mass factorization. We integrate 
numerically over the hard and non-collinear part of phase space and obtain an implicit dependence 
on \og5 s and log<5 c . 

No explicit scale dependence occurs in dimensional regularization after the purely soft diver- 
gences are canceled between the different da s and da sc contributions, and the pre-factors between 
real and the virtual gluon emission diagrams are matched. All poles in da s and all double poles 
in da sc vanish after the soft divergences are removed, and only single collinear poles 1/e remain 
in da sc . Cancellation of these remaining divergences with the soft mass factorization contribu- 
tion renders a finite virtual gluon emission matrix element with an additional explicit dependence 
on log<5 s . This logarithmic dependence cancels against the implicit dependence of the numeri- 
cally integrated hard non-collinear phase space. An explicit dependence on the factorization scale 
log ix p remains in the universal mass factorization terms. Last, as a slight complication, the same 
soft-collinear phase space configurations, which include an explicit log<5 c factor, also lead to an 
implicit log 5 S dependence after the numerical phase space integration. This dependence again 
cancels with the cutoff dependence of the 2 — > 3 phase space, but it makes the numerical analysis 
tedious fl. 

Both cutoff parameters, 5 S and 8 C , must be small in order that the soft and collinear approx- 
imation of the universal terms be valid. Physical observables should not depend on the cutoff 
parameters appreciably. We have checked in detail that the NLO cross section and the distribu- 
tions approach a two dimensional plateau for sufficiently small values of S s and 8 C . To define the 
soft and the collinear regions of phase space consistently and to avoid double counting, it is most 
convenient to require 8 C <C 5 S and then determine the one dimensional plateau for a fixed relation 
between 5 S and 5 C . We show the behavior of the cross section versus the soft and the collinear 
cutoffs in Fig.Elat the LHC energy. The cross section develops a wide plateau for 5 C < 5 S . Unless 
noted otherwise we use 5 C = 1CT 5 and 5 S = 10~ 3 in all numerical analyses. Checking the two 
dimensional plateau and comparing it with results for the total cross section obtained with a one- 
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Figure 2: Cross section dependence on the cutoff parameters S c and 5 S at the LHC energy. Left: The 
soft cutoff S s is fixed to different values. Right: two dimensional logarithmic dependence on both cutoff 
parameters. 

cutoff method Jl^j . we find that the NLO cross section has a remaining uncertainty of 0.1% — 0.5% 
due to the cutoff dependence and the corresponding numerical uncertainty. 



B. Scale Dependence 

Perturbative QCD calculations introduce an unwelcome dependence on the renormalization 
scale hr and the factorization scale //p. One of the major motivations to perform NLO calculations 
is to reduce this scale dependent theoretical uncertainty in predictions of physical observables. As 
the default renormalization scale, we choose Eq. ©, related to the hard scale M. We identify 
the renormalization scales of the strong coupling and the Yukawa coupling. This central renor- 
malization scale choice leadsjo perturbatively stable predictions for cross sections and branching 

22 ] and as functions of y b ll^sll . The situation is different for the 



fractions as functions of a, 



factorization scale. Two reasons point to a central scale considerably smaller than the hard scale 
M, an optimum choice being fi F = M/5. 

First, we can estimate the factorization scale from the kinematics of the process gg — > btH~. 
The bottom quark factorization scale may be defined as the maximum p T b which is included 
in the 2 — s- 2 ^6 — > tH~ process. For the perturbatively calculated bottom quark density we 
can go back to the process gg — ► btH~ and estimate up to which value the cross section 
shows the asymptotic behavior ai tH - ~ l/pr,&- The hadronic phase space or, more specifically, 
the gluon luminosity cuts off the asymptotic behavior near p™^ x ~ M/5. This behavior can be 
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Figure 3: The scale variation of the cross sections at the LHC for mn = 250 GeV and tan f3 = 30. The 
central values are fjP R = M/2 and fi° F = M/5. The four panels show the leading order cross section, the 
complete NLO cross section, the gb induced NLO cross section and the gg induced contribution to the NLO 
cross section. 

understood independently from the form of the matrix element, as long as both incoming partons 
are either gluons or bottom quarks Jl^l fisll . From basic principles it is not clear if one could 
use the (factorizing) internal momentum transfer Qb instead of pr,b llal- The difference in the 
maximum value up to which the asymptotic form holds is p™l x ~ Q b ria,x /2, and the plateau in p T ,b 
is considerably softened [15]. We take this uncertainty into account by varying the factorization 
scale over a generous range. Our argument works because the bottom parton density is calculated 
perturbatively, meaning that its features are well defined and understood perturbatively. 

Second, in the similar LHC process bb — > h, the explicit next-to-next-to-leading-order (NNLO) 
corrections are perturbatively most stable for, and therefore point to, the same small factorization 
scale ifisl . Moreover, the NLO corrections to the similar process bb — > W + H~ are negative for 
\ip = M I*!, indicating possibly a collinear subtraction much too large. 

The variation of the total cross section with the factorization and renormalization scales is under 
control perturbatively if the two scales fx R and fx F are varied independently However, there 
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is a very large shift in the total rate if the two scales are varied together \ir oc \ip and run to very 
small values \i < M/10. This behavior suggests the presence of large contributions proportional 
to log x log up- In Fig. |3] we show the scale dependence of the different contributions to the 
leading order and next-to-leading order cross sections for the process gb — > tH~. The leading 
order curve (upper-left) behaves as one would expect, namely the cross section increases for small 
Hr and for large /if, independently of each other. This behavior arises from the running strong 
coupling a s and the bottom parton distribution. In contrast, the running bottom Yukawa coupling 
is relatively constant for these large scales. The NLO gb initiated curve (lower-left) is also easy to 
understand. The cross section increases with the (bottom quark) factorization scale, but it develops 
a maximum as a function of (jlr at a physical value of the scale, for values of /ip not too large. The 
fraction onlo^/olo * s under control perturbatively over the range of scales, varying at most 20%. 
The corrections are largest (and negative) for small fiR, because a s is largest there. Comparison of 
the NLO and LO results in the OS scheme with those obtained in the MS scheme shows that use 
of the bottom quark pole mass as the Yukawa coupling produces perturbatively less stable rates, 
giving an unacceptably large leading order cross section. The gg — > btH~ contribution (lower- 
right) to the NLO rate has a very different behavior, which dominates the complete NLO rate. The 
gg contribution is regulated through mass factorization, meaning that we compute this process by 
subtracting out the contribution which is included in the collinear bottom quark splitting. For a 
central hr ~ fj° R , the gg channel gives a small positive contribution for small ji F and a small 
negative contribution for large ^f- The latter indicates that a choice of a large scale overestimates 
the logarithmic terms, an overestimate then corrected by the explicit NLO diagram. If we do not 
take into account the difference in size of the gluon and bottom parton luminosities, the gg initiated 
process is suppressed by a factor a s as compared to the gb process. The pattern of correcting 
behavior stays but becomes much steeper if we decrease [ir and thereby increase a s . For a central 
value 11 f ~ l-i'F me 99 initial state corrections are zero, and for much larger or smaller fi F their 
absolute value increases sharply. The gg contribution grows for small /ip and simultaneously 
small jiR. This effect can be rationalized if we take the position that the choice of very small 
Hf and very small /ir corresponds to a region in parameter space where the gb initial state is 
not dominant perturbatively. The result arises because we choose to ignore the presence of large 
collinear logarithms up to p™^ x = iip,b and at the same time we push a s to large values. The K 
factor grows, indicative that we should use the process gg — > btH~ as the leading process. 
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Figure 4: Left: The tH~ inclusive cross section at the LHC as a function of tan (3. The solid curve is the 
NLO result, with the scale variations around the central scale p/fio = 1/4 — 4 indicated by the thinner solid 
curves. Also shown are the leading order result with a running bottom quark Yukawa coupling (dashed 
curve) as well as a pole mass Yukawa coupling (dotted curve). Right: The total cross section at the LHC 
as a function of the charged Higgs boson mass as the solid curve. The NLO and LO cross sections are 
shown as the dashed upper and lower curves, respectively. The dotted curves show the cross sections for 
pp — > tt* + X with a subsequent decay t* — > bH~ in the Breit-Wigner approximation (obw) an d that 
including the complete set of off-shell diagrams (o- ff-sheil)- 

C. Total cross section at the LHC 

The effects of the NLO corrections on the total cross sections for the process gb — > tH~ at the 
LHC are shown in Fig.HJ versus tan/3 and versus the charged Higgs boson mass m H . The solid 
curve shows the NLO cross section, with the scale variations around the central scale /i//io = 
1/4 — 4 indicated by the thin solid lines in the left panel. The total cross section increases for large 
values of tan (3, as expected for the leading behavior cr to t oc tan 2 (3, while the K factor is fairly 
independent of tan j3 \ At leading order we compare the cross section predictions for running 
bottom quark Yukawa coupling and the pole-mass Yukawa coupling. The pole-mass Yukawa 
coupling overestimates the total production rate, while the leading order rate with a running bottom 
quark Yukawa coupling is a more reasonable approximation. With the running bottom quark 
Yukawa coupling, the NLO result is slightly enhanced by a factor of K ~ 1.4 11411 . The scale 
variation suggests a remaining theoretical uncertainty of about 20% on the predicted NLO cross 
section. Both results confirm that the perturbative behavior of the production process gb — > tH~ 
is under control. 
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In the right panel of Fig.|4j the NLO and LO cross sections are indicated by the dashed upper 
and lower curves, respectively. The size of the NLO corrections is essentially independent of the 
Higgs boson mass, a uniform enhancement factor of about 1.4. The process gb — > tH~ itself is 
well defined over the entire range of Higgs boson masses, as long as the hard scale M = m H + m t 
is sufficiently large to motivate the bottom parton picture. 

An interesting region is one in which the Higgs boson mass becomes similar to or smaller than 
the top quark mass, and the decay t — > bH~ is possible. In Fig.|4|we see that the production of a 
top quark pair pp — > ti + X with subsequent (off-shell) decay i — > bH~ becomes the dominant 
process. For m# < m t the production cross section is of order o? s y\ t instead of OL s y\ t , and the 
large gluon luminosity at the LHC is effective in this range of partonic energy. In the following 
we discuss how these two processes can be combined, to obtain a prediction of the cross section 
over the entire range of Higgs boson masses. 

For small Higgs boson masses, below the threshold t — > bH~, the ti production process with a 
subsequent decay of the i dominates the rate for associated production of a charged Higgs boson. 
It is straightforward to combine the ti production process and the exclusive production channel 
pp — > btH~ + X with a tagged final state bottom-quark jet [29]. We compute the process pp — > 
btH~ + X with a finite top-quark width, essentially giving us a Breit-Wigner propagator for the 
intermediate i. Approximate gauge invariance can be achieved in the overall factor scheme [30], 
in which additional terms 0(T t /m t ) are traded for a gauge invariant matrix element. Referring to 
the dotted curves in the right panel of Fig. |4j we see that the process i — > bH~ is approximated 
well by the subset of diagrams with an intermediate i in the Breit-Wigner approximation, as long 
as the Higgs boson mass is below the top quark threshold. Above the top quark threshold, the 
exclusive process pp — > btH~ + X is dominated by the continuum off-shell diagrams, i.e. the 
off-shell extension of the ti production process, including all the diagrams initiated from qq and 
gg. These off- shell diagrams become dominant where the cross section flattens and settles below 
the gb rate, while the Breit-Wigner cross section becomes very small. 

For Higgs boson masses above threshold, a bottom quark jet tag is a heavy price to pay, and it 
is likely a better idea to consider the bottom-inclusive process. The collinear logarithms become 
large, and to obtain the best possible prediction of the rate we must compute the process gb — > 
tH~ [14], keeping in mind that even at threshold the hard scale of the process is M = 350 GeV ^> 
m b . The difference between the bottom-gluon induced rate and the off-shell curve in Fig.|4]shows 
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this enhancement of the rate due to the resummed logarithmic terms. Strictly speaking, we should 
take into account that off-shell production also includes the quark-initiated channels qq — ► btH", 
which do not contribute to the bottom parton density at leading order. However, they contribute 
only about 10% to the total exclusive rate. The matching of the regions of small and large Higgs 
boson mass in Fig. HI indicates that a combination is needed of the ti production process with the 
process gb — > tH~ . There appears to be no region of Higgs boson masses where the off-shell 
production process gg — > btH~ is the appropriate perturbative description. 

In Fig. IU we show the different ti cross sections with a finite bottom quark mass, while we 
neglect the bottom quark mass for the gb channel. The uncertainty induced by this approximation 
is small, however, since we cannot avoid neglecting Y t /m t corrections and rrib/M ~ T t /M. 

At leading order the Breit-Wigner approximation of the pp — ► tt* +X process, with subsequent 
decay of the off-shell t*, may be combined without problems with the process gb — > tH~. We 
may add the independent event samples for any Higgs boson mass value. However, at NLO and 
for Higgs boson masses smaller than the top quark mass there is a potential problem of double 
counting. The ti* production process with a subsequent decay of the t* can be regarded as an 
0(a s ) correction to the gb initiated process, Eq. ©, while it can as well be viewed as (nearly) 
on-shell ti production with a subsequent decay i — > bH~. To avoid double counting, we subtract 
the resonant on-shell part of the ti diagrams from the NLO correction to tH~ production and keep 
it as part of the pp — > ti + X rate. The non-divergent off-shell contribution of the i* propagator is 
counted toward the NLO tH~ rate. The division into on-shell and off-shell contributions, however, 
is well defined only in the narrow width approximation, so we neglect terms of order T t /m t . The 
ambiguity reflects the unsolved problem of how a long-lived intermediate particle is treated in field 
theory. Up to finite width corrections we can, just as at leading order, add the rate for ti production 
with a Breit-Wigner propagator and the properly subtracted NLO gb — > tH~~ rate to obtain a 
prediction for any given Higgs boson mass. In Fig. El we see that addition of the cross sections is 
essentially equivalent to a naive matching procedure. It is perhaps unexpected that the corrections 
from the Breit-Wigner propagator extend to large Higgs boson masses. On the other hand, since 
there is a ~ 20% 1 14] theoretical uncertainty on the NLO cross section for the gb induced process, 
the details of this matching/adding procedure are not important phenomenologically. Instead of 
including the on-shell production with the Breit-Wigner propagator for Higgs boson mass values 
above 300 GeV, we could as well have cut it off at m H ^ m t + 10r t . 
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Below the t — > bH threshold, the Breit-Wigner description is valid, with higher order contri- 
butions included in the cross section for pp —>■ ti + X 113 lL 13211 . Off-shell effects have little impact 



there. Above threshold these off-shell effects have a considerable impact relative to the Breit- 
Wigner description, but the gb — ► tH~ cross section is dominant. After the on-shell contributions 
are subtracted from the NLO rate for gb — > tH~, we can match the two results simply by adding 
them, without any problem of double counting. 

It is useful at this point to compare our final predictions with those published in Ref. [ 19]. For 
m H = 250 GeV and m t = 175 GeV, and for the same value of tan (3 = 30, our predicted K factor 
is 1.4 versus about 1.6 in Ref. J 1911 . However, these numbers should not be compared directly since 
the factorization and the renormalization scales are different in the two calculations. In addition, 
there are small differences in the values of parameters such as the NLO m b , a s and the choices 
of cutoffs. If we use exactly the same parameters, cutoff choices, renormalization schemes, and 
most importantly the same factorization and renormalization scales, our result is 5% larger when 
compared to that of Ref. Jl^l . We attribute this 5% difference to uncertainty in the numerical 
integration. In the two-cutoff phase space method, the final NLO cross section is the difference of 
two large quantities. The numerical uncertainty of either of the these integrations is less than 1%, 
but an uncertainty of 5% can develop in the difference. 

There are two masses in the matrix elements that must be renormalized, m& and m t . The top 
mass m t enters both as an external quark mass and in the Yukawa coupling. The bottom mass m& 
appears only in the Yukawa coupling, since we have set the external bottom quark mass to zero. 
We consistently use the on-shell (OS) scheme for the external top quark mass and we use the MS 
scheme for the top and bottom quark masses in the Yukawa couplings. To understand the effect 
of a different renormalization scheme for top quark mass in the Yukawa coupling, we perform our 
calculation in both schemes. The counter terms for the top quark mass are 

^--^(1^)1, (7) 



m t 4tt e uv 



^ = _^,(fVYr<i +e) (JL + 4), (8) 

m t Aix \ mf J \e uv J 

for the MS and the OS scheme respectively. Changing from one scheme to another can induce 
about 12% difference in the NLO cross section relative to the LO Born cross section. 

We judge it physically more attractive to use the OS scheme for the external top quark mass 
since the top quark mass reconstructed in experiments is the pole mass. Nevertheless, it is perhaps 
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Figure 5: The cross section for the Tevatron (2 TeV) as a function of the charged Higgs boson mass. We 
also show the cross sections for pp —> ti* + X with a subsequent decay t* — ► bH~ in the Breit-Wigner 
approximation and including the complete set of off-shell diagrams. 

best to admit that the difference in scheme choice is tantamount to a difference at the next order in 
a s and thus should be viewed as a theoretical systematic uncertainty at the order of perturbation 
theory in which we are working. 



D. Production at the Tevatron 



The successful matching of the Breit-Wigner approximation and the process gb — > tH~ does 
not apply readily at the Tevatron. In Fig. [5] we observe that the LO process gb — ► tH~ underesti- 
mates the cross section compared to btH~ production 112911 . At the Tevatron the gluon luminosity 
is not dominant in the relevant region of partonic fractional momentum x. The gluon initial state 
contributes only about 10% to the total pp — > btH~ + X rate. Because initial state gluons are the 
dominant source of the bottom partons, we expect the leading order gb rate to be far smaller than 
the true rate. The leading-order 2^3 processes contribute to the NLO gb process, and perturba- 
tion theory for the gb process is not well defined, in the sense that the NLO corrections are large. 
The difference between the leading order gb rate and the off-shell btH~ production rate is slightly 
less than a factor 10 because the gb rate is still enhanced by the resummation of large logarithmic 
terms in the bottom parton picture. We limit ourselves to LHC results for most of the rest of this 
paper. 

The NLO inclusive tH~ rate consistently includes the whole set of quark initiated processes, 
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Eq. ©. Because these quark processes are dominant at Tevatron energies, the K factor for the gb 
process can be as large as 5 for m# ~ 175 GeV. The NLO inclusive rate matches the exclusive rate 
from pp — > t*t + X; i* — > bH~ fairly well, particularly in view of possible remaining differences 
between these two channels. The cross section <r fj_ s h e ii for pp — ► t*t + X is calcualted at leading 
order and evaluated with leading order running couplings and parton densities, while the quark- 
initiated contributions to the NLO inclusive rates cri nc i,NLO are evaluated with NLO quantities. 
Moreover, the collinear divergences in the exclusive rate are regulated by a physical bottom quark 
mass, while the NLO inclusive rate neglects the bottom quark mass and is regulated by mass 
factorization, i.e. by subtraction of the divergent contributions to avoid double counting with the 
NLO evolution of the parton densities. 

Just as for the LHC we see that the Breit-Wigner approximation and the complete off-shell 
matrix element evaluation agree very well up to m H ~ m t . Above threshold the LO gb rate is 
significantly smaller than the complete off-shell rate, but the NLO inclusive rate matches the ex- 
clusive rate well. The visible effect which the Breit-Wigner contribution has on the matched/added 
sum of the cross sections may be unexpected, but we keep in mind the substantial theoretical un- 
certainty on the NLO prediction. The NLO contribution from the qq — > btH~ production process 
is larger than the LO gb induced rate. The formally NLO gb induced rate enters with a much wider 
band of uncertainty than the 20% we quote for the perturbatively well behaved LHC process. This 
wider band covers different schemes for phasing out the Breit-Wigner contribution toward large 
Higgs boson masses. In Fig. 13 we cut off the Breit-Wigner cross sections for Higgs boson masses 
between 200 GeV and 250 GeV, i.e. roughly 20 top quark widths above threshold. 

Although we obtain a predicted cross section at the Tevatron for the entire range of m#, we 
emphasize that the matching of the Breit-Wigner production process and gb fusion works only if 
we take into account the NLO corrections to the gb channel. At the LHC the same matching of 
the two approaches at threshold makes sense even for the leading order gb — ► tH~ rate. Because 
the bottom parton picture is perturbatively stable at the LHC energy, the prediction of the charged 
Higgs boson production rate suffers from smaller theoretical uncertainty. 

IE. KINEMATIC DISTRIBUTIONS 

The bottom parton picture underlies the calculation of the production cross section. As summa- 
rized above, the bottom parton description provides an appropriate way to compute the total cross 
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section for charged Higgs boson production if an appropriate bottom quark factorization scale is 
used. We establish three results in extending our analysis to the kinematic distributions. The nor- 
malized distributions of sufficiently inclusive variables do not change significantly from a LO to 
a NLO treatment of the process gb — ► tH~ . Second, the distributions as well as the total rate do 
not have a strong scale dependence. In particular, checking the limit fi F — > m b , we verify that the 
bottom parton picture does not have much impact on the shape of the kinematic distributions of 
the heavy final-state particles. Finally, we test the approximation of vanishing bottom quark mass 
in the phase space and the matrix elements. 

The first statement is easy to confirm. We show the rapidity, the transverse momentum, and the 
invariant mass distributions for the heavy final-state particles in Fig. |f)J Although the additional 
bottom quark jet in the final state absorbs part of the momentum from the incoming partons, the 
NLO transverse momentum distributions are minimally harder. The extra NLO purely quark- 
initiated production process has a considerably harder p T spectrum, but it contributes only 10% 
to the NLO rate (see also the quark induced contribution in the left panel of Fig. |71). For single 
particle spectra at the LHC, we conclude that the shift in the final state distributions from LO to 
NLO is smaller than the typical scale uncertainty of 20% on the NLO rate. Similar behavior is 



found, for example, in the production of heavy supersymmetric particles at NLO [ 26 



33]. 



A. Zero transverse momentum approximation 

The transverse momentum and the rapidity distributions of the heavy final-state particles are 
depicted in Fig. |6j The choice of the bottom quark factorization scale ixp — M/b has 
to be a part of a consistent bottom parton picture for this class of processes at the LHC 
It remains to be checked whether the collinear approximation for the gluon splitting into a bottom 
parton is appropriate for the distributions of the final state particles. 

We make use of the method described in Sec. Ill Bl to test the bottom parton approach and in 
particular the approximation of negligible transverse momentum of the incoming bottom parton. 
The bottom parton density vanishes if the bottom quark factorization scale approaches the bottom 
quark mass np,b m b- For consistency reasons we use the same factorization scale for all partons, 
but neither the gluon density nor the light quark densities change dramatically as fxp — > rn b ^> 
Aqcd- In this limit the NLO gb induced cross section at the LHC is dominated by the process gg — > 
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Figure 6: The kinematic distributions for the heavy final-state particles at the LHC: transverse momentum 
and rapidity of the charged Higgs boson (left) and the top quark (right) for the process gb — ► tH~. The 
dashed curve shows the leading order distribution with the central choice of scales. The solid and the 
dotted curves represent the NLO results for three choices of the bottom quark and the gluon factorization 
scale fiF = M/5, \ip = 20 GeV and [If = 10 GeV. The bottom row shows the distribution of the 
invariant mass of the top quark and the Higgs boson pair from the LO and NLO calculations, plus a rescaled 
NLO distribution (the last panel). The scaling factor involves the gluon partem densities at the different 
factorization scales. 

btH~. Even though the final-state bottom quark is massless in the calculation, the corresponding 
rate is finite and well defined for any factorization scale fi F > m b . All divergences in the gg — > 
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btH~ process, regulated originally by a bottom quark mass, are absorbed into the definition of the 
NLO parton densities. As described in Sec. Ill Bl when the bottom quark factorization scale takes 
the limit ixpfi —> m b the physical picture shifts from the resummed cross section, including a large 
logarithmic term logpr,b/ m &' to me btH" situation in which the large logarithmic contribution 
could be removed, for example, by a detector cut p™™. We show the difference between the NLO 
distributions with the central factorization scale \i F — M/5 and the small scale limit in Fig.|6j The 
scales we use are ji F = 10, 20 GeV. In principle, we could as well try fi F = 5 GeV, for which the 
calculated bottom quark density is well defined, but the parton densities of gluons and light quarks 
are poorly constrained. We checked that we would then see all the features described below for 
scales down to 10 GeV, except that their effect on the cross sections would be numerically more 
pronounced. 

We see in Fig.|6]that the final-state top quark and Higgs boson momentum distributions become 
somewhat harder when we increase the contributions from the 2^3 matrix elements, going to 
Hf = 20 GeV and /i^ = 10 GeV. The same behavior is seen in the invariant mass of the top quark 
and the Higgs boson pair, m tH , shown in Fig.|6j To explore this feature we present a normalized 
distribution of the m tH invariant mass, in which we rescale the gluon distribution function for the 
central choice /ip = M/5 by a factor P 9 (^f, x PDF ) 2 /P g (M/5, x PDF ) 2 , to estimate the effect of the 
parton densities. Because most of the cross section arises from production near threshold, we can 
approximate xpdf = ram/ E co \\ with E co n = 14 TeV. The physics motivation of this cross check 
is that the gluon densities become slightly harder for smaller scales, and we want to understand 
whether the hardening of the m tH distribution arises from the shift from the bottom parton picture 
to the gluon fusion picture or if it is due to an overall hardening of the gluon parton spectrum. 
In Fig.|6]we observe that scaling the usual NLO distribution with the x dependence of the gluon 
parton density reproduces the hardening of the top quark and Higgs boson spectra. Concluding 
this argument, we find a slight shift in the spectrum at smaller factorization scales, but this shift 
is induced by the shape of the gluon parton density. The two sets of distributions in the bottom 
row of Fig. |6] show that there appears to be no problem with the bottom parton approximation for 
sufficiently inclusive distributions of the heavy final state particles at the LHC. While there is a 
large logarithmic term logpT,&/ m '> present in the btH~ production rate, no additional perturbative 
pitfalls appear in the gb process after the logarithmic terms are resummed. All additional depen- 
dence on the (neglected) transverse momentum of the final-state bottom jet is effectively power 
suppressed. 
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Figure 7: The kinematic distributions for the final state at the LHC: the top quark transverse momentum 
and the invariant mass of the top quark and the Higgs boson pair. We compare the NLO result for the 
process gb — ► tH~ with the process pp — ► btH~ + X. Apart from the complete set of diagrams with the 
physical bottom quark mass, we show the purely qq induced process and the complete set of diagrams with 
a mathematical cutoff instead of the physical bottom quark mass m& — ► 0.46 GeV. 

B. Zero bottom quark mass approximation 

We have compared the gb induced process with the process gg — ► btH~ without a finite bottom 
quark mass in the phase space or in the matrix element. For the kinematic distributions of the 
bottom quarks this approximation is not obviously good. In our NLO approach the divergences in 
the p T b spectrum are compensated by a negative infinity at pr,b = 0, i.e. in the 2 — > 2 kinematic 
limit. This distribution is not physical, and all-orders soft gluon resummation should be taken 
into account [34] to obtain a physical spectrum with a peak at some small value of pr,b- With the 
bottom quark mass as a regulator, the px,b spectrum peaks near m b [ 14]. However, when the gb 
process is used, we are implicitly not interested in observing the final bottom-quark jet and in its 
distributions; rather, we are interested in the distributions of the heavy final state particles. In Fig.|7] 
we show the normalized transverse momentum distribution of the top quark for the gb process at 
NLO, and for the 2 — > 3 process, with two different cutoffs: one with the physical bottom quark 
mass and the other with a smaller mathematical cutoff (we use 1/10 of the bottom quark pole 
mass). We observe that the gb calculation agrees with the 2^3 matrix element approach with 
the physical bottom quark mass. The curve with a smaller cutoff instead of the bottom quark mass 
agrees perhaps too well with the NLO process in which the bottom quark mass is neglected. The 
distribution in the invariant mass of the tH~ final state confirms this level of agreement. The 
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dependence on the bottom quark mass seems to be power suppressed, and the approximation of 
zero bottom quark mass is justified. 

There are limitations of our argument for charged Higgs boson production at the Tevatron. The 
2^3 rate with a finite bottom quark mass, induced by incoming quarks, shows considerably 
harder momenta of final-state top quarks. While the total cross section is predicted correctly in the 
bottom parton picture, this effect might mean that the kinematic distributions require more careful 
study at Tevatron energies. The effect is not related to the approximation of zero bottom quark 
mass. Instead, it probes the perturbative link between the contribution of gluon and quark initiated 
diagrams in the 2 — > 3 production processes and the bottom parton description. 

Finally, we remark that the total cross section and the inclusive distributions of the final state 
particles are correctly predicted in the bottom parton picture. Neither the small transverse momen- 
tum approximation nor the small bottom quark mass approximation in the bottom parton picture 
has a visible effect on the transverse momentum, the rapidity, and the invariant mass distribu- 
tions of the final state top quark and Higgs boson. Shifts induced by the massless bottom parton 
approximation are washed out once detector resolution is taken into account. For the processes 
under consideration, the contribution from the qq initial states is less important at the LHC than at 
the Tevatron. The light-quark induced subprocesses show harder transverse momentum spectra, a 
difference that should be considered for predictions at Tevatron energies. 

C. Final state correlations 

The fully differential nature of the two-cutoff method enables us to place a kinematic cut on 
one final state particle and to study the distribution in the momentum of other particles in the 
final state. It allows us to examine any correlation observable among the final state particles 
which does not spoil the cancellation of soft and (initial state) collinear divergences. For reasons 
of acceptance and/or to improve signal purity with respect to backgrounds, it may be helpful to 
make final-state cuts that act in similar fashion to a cut on the transverse momentum of the top- 
quark. If the charged Higgs boson decays to a r-lepton jet, a simple transverse momentum cut 
on a lepton from the top quark decay could be an example. We examine in this section how 
such a selection may affect the expected transverse momentum distribution of the Higgs boson. 
The Higgs boson transverse momentum distribution is crucial because it determines the boost of 
the Higgs boson decay products. If these products include bottom-quark or r-lepton jets, their 
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Figure 8: LO and NLO p^-distribution of the Higgs boson, with cuts px.t > 50 GeV or pr t t > 100 GeV 
on the transverse momentum of the top quark. The Higgs boson mass is uih = 250 GeV, and tan f3 = 30. 

transverse momentum distributions determine the tagging efficiencies. 

The topology of the final state and the correlation in momentum between the top-quark and 
the Higgs boson change once NLO corrections are included. Instead of a back-to-back pair of 
heavy states at leading order, the NLO topology includes three final-state particles sharing the 
total transverse momentum. In Fig. HI we show the pr,H distributions after the cuts px,t > 50 GeV 
or px,t > 100 GeV. At LO, a cut on px,t indeed eliminates values of pt,h below this cut. The 
LO distribution da/dpT,H in the range pt,h > 100 GeV is identical for the two cases shown: 
p™j _ 5Q Q e y and p min _ 1Q q Ge y At nlo, the figure shows that the effects of cuts on the 
momentum of one final state particle extend over a significant range in the momentum of another 
final state particle. Owing to the presence of the final state jet, the NLO transverse momentum 
distribution of the Higgs boson extends all the way to zero, as shown in Fig.[U The impact of the 
PT,t cut on the NLO distributions is evident well above pt,h = Prt> an( ^ me distributions do not 
coincide until pt,h > 200 GeV. The ratio of the NLO and LO distributions cannot be represented 
by a simple correction (or "K") factor. This factor would be infinite for values of pt,h less than 
the value of the cut on p T)t , less than unity in a small interval where p t>H is just above the cut, and 
uniformly greater than unity for p TjH > 1.5p™ n . 

In contrast to the differences seen in Fig.[H the effect of a pr,t cut on the Higgs boson rapidity 
distribution is trivial: the NLO rapidity distributions of the Higgs boson and the top-quark are 
reduced by an approximately uniform factor. 

The distribution in p T>H is well defined at NLO. The soft and the initial-state collinear di- 
vergences appear with Born-type kinematics, just like the explicit infrared poles. However, we 
expect strong cancellations between large NLO negative virtual contributions and positive NLO 
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real emission contributions at the LO threshold. Numerical problems can arise once these cancel- 
lations yield a cross section at threshold that is more than an order of magnitude smaller than the 
contributing parts. Motivated by evidence of instability in the threshold region p T:H = p™\ if we 
use the values of 5 S and 5 C chosen earlier, we select cutoff values S s = 1CT 2 and 5 C = 2 x 10~ 4 
to obtain the results shown in Fig. [5] These relatively large values of the cutoffs remain in the 
safe region for the total cross section, as shown in Fig. |2j Moreover, we use fairly wide bins, to 
be less sensitive to numerical cancellations in the threshold region. The behavior of the curves 
in Fig. |E] near threshold indicates a some remaining uncertainty, so the p TH distributions should 
be taken with reservation at threshold. They should be reliable one bin or more from threshold in 
each direction. 

Another distribution of interest is the vector sum of the transverse momenta of the top quark and 
the Higgs boson (£pr = | (pr,t + Pt,h)\) shown in the left panel of Fig.[9j The LO distribution is a 
single peak at Hp? = since the top quark and the Higgs boson balance in transverse momentum 
at this order. At NLO, the observable Ep T is the transverse momentum of the third jet in the final 
state. The shape of the distribution depends only mildly on the Higgs boson mass. The NLO 
distribution in Epr shows a marked divergence as T,px — > 0. This divergence reflects a limitation 
of our fixed-order calculation and points to the eventual need for all-orders resummation of the 
effects of soft-gluon radiation, as discussed and implemented for other processes [34]. 

The right panel of Fig. |9] displays the two-dimensional correlation between Hp? and 
A(j) = \(fi t — (j) H \ for the sum of all 2 — > 3 contributions to the cross section. Here is the 
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Figure 9: Left: Normalized distributions in T,px (defined in the text) for the NLO process gb — > tH , for 
niH = 250 GeV and 500 GeV. Right: A two-dimensional plot of the correlation between Sp^ and A<f), 
for m H = 250 GeV and tan/3 = 30. 
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Table I: Supersymmetric corrections to the production cross section gb — > tH~ from non-resummed and 
resummed A^, corrections, Eq. ©, and from the explicit remaining supersymmetric loop diagrams. The 
supersymmetric parameter points are chosen according to the benchmarks in Ref. [38]. All masses are 
given in units of GeV. The percentage changes are defined with respect to the purely gluonic NLO rates. 

difference between the azimuthal angles of the Higgs boson and the top quark. At LO it is fixed to 
Acf) = 7T, but for the 2^3 processes the distribution extends to all A0, again a manifestation of 
the fact that an additional parton in the final state contributes to the transverse momentum balance. 
It is interesting to see that the Ep T distribution for A</> ^ ir, i.e. for a topology different from LO, 
develops a maximum for Epr ~ 50 GeV and drops to zero for T,px — > 0. This feature and the 
corresponding behavior of A0 for T.px ^ reflect the fact that the (resummable) divergence is 
limited to the LO topology. 

IV. SUPERSYMMETRIC CONTRIBUTIONS 

Supersymmetric diagrams contribute to the production rate for gb — > tH~ at the same level 
as the NLO QCD contributions (a^y^), These diagrams are virtual gluon exchange diagrams, 
where the gluons and quarks are replaced by their supersymmetric partners, gluinos and squarks. 
A feature of squarks is the mixing between the supersymmetric partners of the left-handed and the 
right-handed quarks. The 2x2 bottom-squark mass matrix has an off-diagonal entry m b (A b — 
\i tan (3), and the top-squark mass matrix has an entry m t (A t — ji/ tan 0). Here fi is the Higgsino 
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mass parameter which links the two Higgs doublets in the Lagrangian, and A4 is the trilinear 
squark-squark-Higgs boson coupling parameter. The off-diagonal entry in the matrix element is 
proportional to the mass of the standard model partner, and it is usually neglected for the first and 
second generations. For large values of tan (3 bottom squark mixing can become larger than top 
squark mixing. 

In the limit of large tan (3, the leading supersymmetric contributions are not loop diagrams, but 
renormalization terms 13511 . This point becomes evident if we compute the corrections to the btH~ 
vertex in the limit of vanishing bottom quark mass, but finite bottom quark Yukawa coupling. This 
approach is justified from a formal point of view because the connection between the mass and 
the Yukawa coupling is a property of electroweak symmetry-breaking and not protected once the 
symmetry is broken. Disassociation of the mass and the Yukawa coupling becomes apparent in 
a type-II two-Higgs-doublet model like the MSSM. Because the bottom quark Yukawa coupling 
always appears as m;, tan (3, the relationship between the mass and the Yukawa coupling is not 
fixed. If we compute the renormalization of the tbH~ vertex with a zero external bottom quark 
mass, bottom squark mixing diagrams lead to contributions which look like mass renormalization 
terms, i.e. terms which create a finite bottom quark mass in the external leg. However, this inter- 
pretation cannot be correct, since mass renormalization has to be multiplicative. Instead, we see 
that these renormalization factors describe a misalignment of the bottom quark Yukawa coupling 
and the bottom quark mass, which appears for zero bottom quark mass as well as for finite values. 
In complete analogy to a mass renormalization, these coupling renormalization diagrams can be 
resummed to all orders and lead to a correction ihfll : 

nib tan (3 m b tan (3 1 
v v 1 + Af, 

Ab= Sm lT^ £ ° F m§ ^ [ B (°' m b^ m §) -B{0,m- b>1 ,mg)] 



-^-C F m g (-A b + /itan/3) I(m~ bl ,m~ b2 ,mg) 



I (a, b, c) 



a 2 b 2 c 2 ' 

a 2 b 2 log — + b 2 c 2 log — + cV log — 

b z c z cr 



(9) 



(a 2 - b 2 )(b 2 - c 2 )(c 2 - a 2 

The function B(p 2 , mi, m 2 ) is the usual scalar two-point function; CV = 4/3 is the Casimir factor 
in the fundamental representation of SU (3). There are similar additional terms proportional to the 
strong coupling or to the top quark Yukawa coupling, but Eq. © is the leading contribution for 
large tan (3. 
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Figure 10: The supersymmetric contributions to the NLO cross section for gb — > tH~, with the dotted curve 
showing the resummed A& contribution and the solid curve showing the combination of the resummed A5 
and non-At contributions. The supersymmetric parameters are described in Table|l] The central value of \i 
in SPS4 is noted in the plot; the central value of tan (3 in SPS5 is at the lower end of the plot, tan (3 = 5. 
All parameters and masses except for fi and tan (3 are kept constant. The curves denoted o"2hdm show NLO 
QCD cross sections without SUSY contributions. 

Since these A b corrections are the leading tan /3-enhanced supersymmetric contributions to the 
production cross section, and since the charged Higgs boson search is most promising at tan (3 > 
15, we might speculate that these corrections to the gb — > tH~ production rate are sufficient. 
Equation © shows that the shift in the Yukawa coupling can have large effects J37I1 provided that 
fx tan (3 is large (preferably negative), the gluino mass is large, and bottom-squark masses are not 
too large. In this limit the percentage corrections are approximately A;, ~ fit&nfl/rrig. The non- 
A b -type supersymmetric corrections are negligible compared, for example, to the remaining NLO 
scale variation. 

The supersymmetric contributions to the charged Higgs boson production cross sections are 



shown in Table[T]for the "Snowmass points and slopes" (SPS) 113 811 parameters. The contributions 
are split into the A& corrections, as defined in Eq. ©, and the remaining supersymmetric diagrams. 
We present the A& contributions in the NLO version 1 — 2A& as well as after resummation, 1/(1 + 
A5) 2 . The ratios of the SUSY corrections to the NLO QCD cross sections are given in the last 
three columns. The negative sign of the Ab contributions is fixed by the sign of fi. The sign of 
is a free parameter, linked to SUSY contributions to the transition rate b — ► sj. In this process 
the measured rate is consistent with the standard model prediction. There are additional charged 
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Higgs boson and chargino induced contributions in the MSSM. For [l > they enter with opposite 
signs and therefore tend to cancel numerically, while for \i < 0, in particular in the mSUGRA 
supersymmetry breaking scheme, the parameter space is closely constrained. Therefore, all SPS 
points are chosen with positive sign of /x. 

In Table U we observe that the A b corrections are dominant for all points with tan/3 > 15, 
particularly for the two points with tan (3 = 30, 50. The leading contribution for large tan (3 is 
described correctly by the A b corrections. At maximum, all supersymmetric corrections are of the 
order of the remaining scale variation and our estimate of the theoretical uncertainty of 20%, as 
long as tan (3 < 30. This modest correction is not necessarily true for the entire supersymmetric 
parameter space, and the A b corrections can be much larger [37]. However, the small correction 
reflects the ansatz used in supersymmetry breaking. None of the scenarios in Tablets designed to 
produce a large splitting in the supersymmetric mass parameters at the weak scale or a large value 
of \fj,\, which would favor large Afe-type corrections. In general, large values of |/i| are a challenge 
in high-scale motivated models. For fine-tuning reasons these models usually produce of order 
the weak scale, to avoid large cancellations of different renormalization group contributions to the 
value of m z . Even in the focus-point [13911 inspired SPS2 the contribution to weak-scale param- 
eters that are proportional to m Q cancels in itself, decoupling the value of m from the leading 
renormalization group running. All other parameters remain at typical weak- scale values. In all 
three SUSY breaking scenarios considered, the large gluino mass, linked to the relative dominance 
of the corresponding beta function (3 3 , is the reason the A b correction is not negligible. 

In Fig.[K)]we show that the A b contributions can become large once we depart from the unifi- 
cation scenarios. Starting from the mSUGRA motivated points SPS4 and SPS5, listed in TableHl 
we vary /j, and tan (3, leaving all other masses and parameters invariant. As expected, the non- 
resummed A b corrections become arbitrarily large for large values of and the resummed A b 
correction can become arbitrarily large for some negative values of /i, both limited only pertur- 
batively and ultimately by unitarity of the enhanced bottom quark Yukawa coupling. The sign 
of the A b correction is fixed by the sign of fi, and the remaining SUSY contributions are small. 
For comparison, we also show the NLO cross sections without SUSY contributions included, the 
typical case for a two-Higgs-doublet model (2HDM). These curves are labeled ct 2 hdm- The NLO 
QCD cross section does not depend on the SUSY parameter /i, but it increases as tan (3 gets large, 
as expected. In Fig.[TT]we show the same effect, starting from the scenarios A and C in Ref. [13711 . 
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Figure 11: The supersymmetric contributions to the NLO prediction for the cross section gb — > tH~. 
The supersymmetric parameter points are the scenarios A and C, picked from Ref. I37L All parameters 
and masses except for \i are kept constant. The gluino mass in both scenarios is 1 TeV; the lighter top- 
squark and bottom-squark masses are 0.5 TeV for scenario A (LHS) and 1 TeV for scenario C (RHS). The 
top-squark mass difference is 100 GeV and the bottom-squark mass difference 150 GeV. In contrast to 
Ref. [37] we resum only the SUSY-QCD corrections. The curves denoted 02HDM show NLO QCD cross 
sections without SUSY contributions. 

Because the value of tan /3 = 50 is large, the non-Aj, corrections are completely negligible, while 
the A b corrections can become arbitrarily large. The difference in the size of the corrections in 
the two panels of Fig. [TT]can be understood from Eq. © in the limit a ~ b for the bottom squark 
masses and either c a, b or c ~ a, b. In both cases the A b corrections are suppressed by the 
heaviest mass in the system, i.e. the gluino mass, but the pre-factor is larger if all masses involved 
are of the same order. 



V. CONCLUSIONS 

We evaluate the inclusive and differential cross sections for the associated production of a top 
quark along with a charged Higgs boson at Tevatron and LHC energies to next-to-leading order in 
QCD and in supersymmetric QCD. 

Using the two-cutoff scheme to treat the soft and collinear singularities, we find stable results 
for total and differential cross sections over large ranges of the cutoff parameters as well as of the 
factorization and renormalization scales. While the QCD corrections to the total rate are sizable at 
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theLHC,AT~ 1.4 [..14], the shifts in the normalized kinematic distributions of the heavy final state 
top quark and Higgs boson are negligible. The scale dependence gives us a reasonable estimate of 
about 20% on the remaining theoretical uncertainty. 

In the regime where m H < m f , we compute the NLO cross section by subtracting the interme- 
diate on-shell divergences in the narrow width approximation. This procedure allows us to match 
the NLO cross section for the process gb — > tH~ with the contributions from gg — ► ti* with a 
subsequent decay t* — > bH~, simply by adding the rates. This method yields a prediction for the 
cross section for associated charged Higgs boson production over the entire range of Higgs boson. 

At the Tevatron, charged Higgs boson production is likely to be observed only for small masses 
of the Higgs boson. In this regime we show the validity of the Breit-Wigner approximation in the 
ti production process. We can add the off-shell production rate at NLO. 

Examining the NLO momentum distributions for inclusive charged Higgs boson production, 
we show the validity of the bottom parton description beyond the total rate. Neither the collinear 
phase space approximation nor the approximation of zero bottom quark mass has a visible impact 
on the kinematic distributions of the heavy final state particles. 

The fully differential nature of the two-cutoff method enables us to place a kinematic cut on 
one final state particle and study the distribution in momentum of the other particles. It allows us 
also to examine momentum correlations among the final state particles. 

We explore the effects of virtual supersymmetric particles in NLO loop diagrams and find that 
the universal A& corrections to the Yukawa coupling can be sizable. In the two Higgs doublet 
model, the remaining explicit loop contributions to the NLO rate are below the level of the scale 
uncertainty. 
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